Quantum walk on circles in phase space 
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We propose a variation of the quantum walk on a circle in phase space by conjoining the Hadamard 
coin flip with simultaneous displacement of the walker's location in phase space and show that this 
generalization is a proper quantum walk albeit over multiple concentric circles in phase instead 
of just over one circle. We motivate the conjoining of Hadamard and displacement operations by 
showing that the Jaynes-Cummings model for coin+walker approximately yields this description in 
the dispersive limit. The quantum walk signature is evident in the phase distribution of the walker 
provided that appropriate pulse durations are applied for each coin flip. 



I. INTRODUCTION 

The quantum walk (QW) is one of the most impor- 
tant developments in theoretical quantum information 
science, both as an intriguing generalization of the ubiq- 
uitous random walk (RW) in physics QL Q and for expo- 
nential algorithmic speed-ups [3, 0, ■ The quincunx, 
or Galton Board Q, was developed to exhibit the fea- 
tures of random walks in experiments and more recently 
an optical quincunx that simulates a 'wave walk' Q was 
demonstrated. 

For the quantum quincunx, an appealing strategy for 
experimental realization arises in the context of a QW 
over a circle in phase space @ , which arises naturally for 
a simple harmonic oscillator. Points in phase space cor- 
respond to the oscillator position- momentum pair (x,p), 
which we henceforth refer to as the phase space 'loca- 
tion', and energy-conserving evolution of the oscillator 
guarantees that E = (x 2 +p 2 )/2 (for the oscillator of 
unit mass and unit frequency) is a conserved quantify, 
thereby constraining the phase space trajectory to circle 
in phase space centered at the origin (0, 0). 

The discrete walk on the circle corresponds to phase 
jumps A6* = 6*2 — 6*1 for 



0; 



tan 1 pi/x i , 



(1.1) 



which is well-defined provided that x ^ ^ p. The 
discrete random walk on the circle, corresponding phase 
jumps ±A6>, with A8 of fixed size and the sign ± chosen 
randomly, has been used to provide a clear explanation 
of phase diffusion of the laser field [l(|. More recently 
the random walk on the circle in phase space has been 
generalized to the QW on a circle in phase space: in the 
quantum case the walker's location as a point in phase 
space is replaced by a localized wave function centered at 
a location (x,p), and the random flip of sign ± is replaced 
by a quantum coin given by a qubit, which is flipped 
by a Hadamard operation and then entangled with the 
oscillator by free evolution. An example of a localized 
wave function is the coherent state 



| a) = D(a)\0) 



(1.2) 



with |0) the ground state of the simple harmonic oscilla- 
tor and 



D(a) — exp(aa^ — a* a) 



(1.3) 



the unitary displacement operator [llj . and a = (x + 
ip)/V2 for localization at (x,p). The full quantum walk 
on a circle in phase space is described in detail in Sec. [Til 

In these quantum-walk-on-the-circle schemes, the coin 
qubit is directly controlled; however, in the context of 
cavity quantum electrodynamics, the coin qubit is an 
atom within a high- finesse resonator. The high- finesse 
nature of the cavity mitigates against direct control of the 
atom, and the only viable option is to drive the coin qubit 
indirectly by driving the cavity, which in turn drives the 
atom. Here we show that this indirect coin flip indeed 
suffices to create a quantum quincunx, but the quantum 
walk is no longer confined to a circle in phase space but 
rather undergoes a quantum walk that hops between dif- 
ferent circles in phase space. In Sec. IIIII we explain this 
revised QW involving simultaneously driving of both the 
oscillator and the coin qubit. 

In Sec. IIVI we consider the Jaynes-Cummings (JC) 
model [l2j as the underpinning of the results in Sec. IIIII 
Whereas Sec. lIIII presents a generalized Hadamard trans- 
formation of the coin that involves driving both the os- 
cillator and the coin, in this section we use the JC model, 
which is one of the most important models in quantum 
optics to describe cavity quantum electrodynamic sys- 
tems, and was originally used to explain the maser (and 
hence the laser) [121 ] . In the dispersive limit, and with ju- 
dicious timing to achieve the right phase steps, we show 
in this section that a QW on circles in phase space can 
be well approximated by JC dynamics. The results are 
summarized in Sec. [V] 



II. BACKGROUND 

The random walk on the circle in phase space, used to 
describe laser diffusion [lofl . comprises two coupled sys- 
tems: the walker, who is physically a simple harmonic os- 
cillator, and the unbiased two-sided coin, which is math- 
ematically an unbiased random bit. The joint system of 
the coin+walker has a state space C 1 (K) x {0, 1}. That is, 
the walker's state corresponds to distributions in £ 1 (K.), 
and the coin can have either value <r € {0, 1}. Evolution 
consists of alternating coin flips, which generates or 1 
randomly with equal probability, and then the walker's 
distribution in phase space is rotated by an angle ±A9 
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with the sign ± given by (— 1)^. 

In quantizing the QW, the walker's distribution is re- 
placed by a state p g B(H W ) for B(H W ) the Banach space 
of bounded operators on H w = £ 2 (R). The coin is re- 
placed by a qubit with Hilbert space H c — PC 2 , namely 
the projective space of two-component complex vectors. 
The joint coin+ walker space Tl c <S> Ti. w is spanned by a 
basis set comprising tensor products of Fock states 



it" 



(2.1) 



for |0) the oscillator ground state in Eq. (|2.1[) . and |0) 
and |1) the two coin basis states. The Fock states are 
also known as number states, and Fock state \n) is an 
eigenstate of the number operator n = a* a with eigen- 
value neN. 

The QW is effected as an alternating sequence of two 
operations, namely the Hadamard transformation on the 
coin 

ff = |+)(0| + |-)(l|, |±> = (|0>±|1»/V2, (2.2) 
and the free evolution 



F(A0) = exp(ina z Ad) 



(2.3) 



between coin flips. The free evolution effects a condi- 
tional rotation of the walker's state by an angle ±A8 
which is chosen given an initial walker state \a) [Hj]: 



1 2tt 
< A9< j=. 



(2.4) 



In fact the evolution of the walker can be entangled with 
the coin state by this evolution, and this entanglement 
between the coin and walker degrees of freedom under- 
pins the dramatic differences between the classical ran- 
dom walk vs the QW. The resultant evolution is achieved 
by repeated application of the QW unitary operator 



U = F(H®1) 



(2.5) 



after N discrete time steps, the state of the coin+walker 
evolves according to the evolution operator U N . 

As we shall see, the QW signature will be evident in 
the phase distribution of the walker's state [HI, 0]. The 
phase distribution for the walker's reduced state p w , ob- 
tained by tracing out the joint coin+walker state over 
the coin's degree of freedom, is 



P{4>) = lim m(0|Pw|<£)m 

M— >oo 



as constructed from phase states [1 

M-l 

i 

<M ~ 



M-l 



II 10 



n) 



(2.6) 



(2.7) 



n=0 



Phase states are thus dual to the Fock states in the sense 
that 



H0) 



M 



/VM 



(2. 



if n < M, and the overlap is zero otherwise, and, for 
4> m = 2rrm/M, 

span{|0 m ) ;m = 0, 1, . . . , M - 1} 

=span{|n) ;n = 0, 1, . . . , M - 1} (2.9) 

with {|</> m )} an orthonormal basis of the subspace. For 
arbitrary phase states \<j>) and \<f> + 5), their overlap is 
given by 



M-i 



(2.10) 



M 



for 



U M -i{5) 



sin(M5/2) 
sin((5/2) 



(2-11) 



the Chebyshev polynomial of the second kind [lq . 

In the coin+walker basis, with phase states as the 
walker basis states, the free evolution operator acts ac- 
cording to 



Fk, 



k,0+(-l) ; A0), ? G {0,1} 



(2.12) 



so the phase states form a natural representation for 
studying this evolution. Furthermore the signature of 
both the random walk and the QW, and their differences, 
is in the phase distribution (|2.6|) of the reduced walker 
state p w . 

The dispersion of the phase distribution is especially 
important. As moments are not particularly useful for 
distributions over compact domains, other strategies are 
needed. For the phase distribution over the domain 
[0, 2tt), Holevo's version of standard deviation [I?) is par- 
ticularly useful as it reduces to the ordinary standard 
deviation for small spreads and is sensible when the dis- 
persion is large over the domain [X 81 ] - Holevo's standard 
deviation is 



2tt 



x /\{j*)\-*-l,(e i +) = j ctyP(0e* (2.13) 

with respect to any phase distribution P{<j>) (|2.6[) . 

The Holevo standard deviation has been shown to 
evolve according to cth oc t for the QW, whereas cth oc \fi 
for the RW, at least for short times where the pha se dis- 
tribution has support over less than the circle [13| . This 
quadratic speed-up of phase spreading in a unitary evo- 
lution is a hallmark of the QW on the circle. We will use 
this quadratic speed-up as the indication of QW in the 
system. 

Our focus is on phase spreading as a signature for 
the QW, but phase is not directly measurable. However 
phase can be inferred from homodyne or from optical 
homodyne tomography measurements (l3j . 
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III. QUANTUM WALKS ON CIRCLES 



In previous schemes, implementations of the QW on 
a circle have been proposed for ion traps Q or cavity 

qed mm 

and each scheme relies on direct driving of 
the coin (i.e. directly flipping the coin without modifying 
the cavity field). In realistic systems this may not be 
possible, and instead the simple harmonic oscillator will 
be driven, which then drives the coin via the oscillator- 
coin coupling. 



A. Generalized Hadamard Transformation 

In this section we treat this strategy of indirectly driv- 
ing the coin by generalizing the Hadamard transforma- 
tion to 



H i ► exp ji^ [a x + Ai] j = H ® £>, 



(3.1) 



for D(a) the unitary displacement operator (|1.3[) with 
a i— ► iX/y/2. Thus A is the kick the walker receives during 
the Hadamard pulse. In the next section we will derive an 
approximation to the unitary operator (|3.1[) by beginning 
with the JC model Hamiltonian. 

The generalized Hadamard transformation (|3.ip nicely 
factorizes into a Hadamard transformation and a dis- 
placement operation. The Hadamard tranformation ef- 
fects the desired coin flip, but the displacement opera- 
tor simultaneously moves the walker to another circle in 
phase space. As A in Eq. (|3.ip is real, the kick is a dis- 
placement in x. The nature of QW on circles in phase 
space is made clear in Fig. [T] 

In this geometric representation, the coin flip 
Hadamard operation is accompanied by a concomitant 
displacement that shifts the walker's distribution (the 
large black dot in Fig. [T]) from one circle of radius rij 
to another circle of radius ny . To understand the effect 
of hopping to different circles of phase space, let us con- 
sider a coin+walkcr state initially in the state |0, a) with 
a = ^^pj which corresponds to the coin in the state 
and the walker localized at {x,p) in phase space. 




Phase spread 



with coordinates x and p, 
Wj for three different values 



FIG. 1: Phase space diagram, 
depicting circles of fixed radius v 
of rij corresponding to three circles indexed by j. The large 
black dot on the innermost circle represents the distribution 
of the walker. The arrows extending from the large black 
dot represent a fixed phase jump ±A#, which is clockwise or 
counterclockwise depending on the sign of AO. The dashed 
line shows a jump size of AO and the corresponding arcs of the 
circles subtended by this angle. The geometric meaning of the 
displacement operator parameter a is shown as the modulus 
being the radius of the j th circle; not shown is that arga is 
the phase of the center of the large black dot. 



state (|3.3p . which yields the resultant state 



U 



0, 



x + ip 



y/2 I a/2 



x + X + ip ^ Ae 
' V2 



, x + X + ip _■ 



AO 



(3.4) 



B. The First Step 

The first step corresponds to the application of the 
unitary operator 



U = F(H®D). 



(3.2) 



First the generalized Hadamard transformation H 
D(iA/v2) is applied: 



x + lp 

~7f 



x + X + ip 
x/2 



(3.3) 



This generalized Hadamard operator is then followed 
by the unitary conditional phase operator F on the 



Eq. (|3.4|) has three important features. One is that the 
resultant state is an entanglement of a coherent state 
with a qubit of the type that is observed in microwave 
cavity quantum electrodynamics experiments The 
second important point is that each of the two walker 
states \(x + X + ip)/v / 2e ±lA6 '} are localized on the same 
circle in phase space, and third the rotation of the coher- 
ent state by angle AO is independent of which circle the 
walker is on. 

Thus, although the walker is forced to hop between 
circles during the application of each Hadamard trans- 
formation l|3.ip . we will show that the QW survives this 
generalized action. 
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The state after N steps is |$(iV)) = U N |$). The phase 
FIG. 2: The phase distribution P(<j>; N, A, M) for the walker's distribution for the walk after N steps is P((j>) in Eq. IpTf]) 
location after A = 4 steps of the QW over the different circles f or Pw = Tr c (|$(A)) ($(A)|). The first 3 steps for the 
in phase space with initial state (|0, a) + i |1, a))/y/2, A<9 = walk and corresponding phase distributions are discussed 
0.35 and (a) A = 0, (b) A = 0.2, (c) A = 0.3 and (d) A = 0.4. in thc Appcndix 



The state after N steps is 

2 N-1 

\$(N,a,\,A6)) = £ \pi(N) |0, oii(N, a, A, A9)} + Qi(N) |1, /3»(iV, a, A, A#)) 



(3.6) 



The coin+walker state (|3.6[) adopts a simple form: it is an entanglement between orthogonal coin qubit states with 
superpositions of coherent states. The weights p%,Pj and coherent state amplitudes ev^, /3j are determined by recursion 
relations presented in thc Appendix. After tracing out the coin state, 

Pw(N, a,\,A6)=J2 [Pi WP* (AO It* (A, a, A, AO)) (a, (A, a, A, A6)\ + Qi {N)q* (A) |^ (A, a, A, A6>)) (A, a, A, A0) | 

(3.7) 

is obtained. The phase distribution for the state after N steps is thus 

P(0;A,a,A,A0) = lim V L(A)p*(A) A/ (0|a l (A, a, A, A0)) x (a 3 -(JV„ a, A, Afl)|0) M 

M— *oo — ' L 



i -J 



+ q l (N)q*{N) M (<f ) \(3 l (N,a,X,A9)) x ((3j(N, a, A, A6*)|0)a/ 
where the overlap of thc phase state with the coherent state given by 



M-l , ia -ich\ n 



(3.8) 



(3.9) 



which is a function of both A and A6. 



D. The Spread in Phase 

In order to observe a QW, the choice of parameters is 
critical. Therefore, we study how choices of AO and A can 
affect the quality of the phase distribution for revealing 
a signature of a QW. We expect that the choice of AO 
controls the rate of spreading of the phase distribution 
because AO corresponds to the size of the walker's step. 
On the other hand, A is responsible for breaking the sym- 



metry of P(<j>) around = 0. We can sec these effects in 
Fig. [21 Specifically, we observe that for increasing A, the 
overall distribution becomes more skewed towards posi- 
tive (f> and individual peaks can become narrower. Thc 
skewing is due to the increasing contribution from |/3j), 
which can be higher in amplitude than the |aj) terms, 
hence the concomitant narrowing of some peaks. 

The spread of the phase distribution provides an im- 
portant signature of the QW, and we use the Holevo 



5 




the random walk and the quadratic enhancement for the 
QW. Therefore, the QW behavior is clearly present de- 
spite having generalized the Hadamard transformation to 
Eq. (|3.1[) and used a Holevo standard deviation for phase 
as a quantifier. Fig. [3] thus makes it clear that the QWs 
over different circles in phase space are actual QWs. 



FIG. 3: The Holevo standard deviation oh of the phase distri- 
bution for the random and quantum walks, for a = 3, A = 0.4, 
and A8 = 0.35, as a function of the number of steps N pre- 
sented as (a)aH vs N and as (b) logon vs logN for the clas- 
sical random walk (solid line) and the QW (dots) . 

standard deviation ch 13f> to quantify this spread. 
The graphs of cth vs logiV and its log- log version in 
Fig. [3] clearly reveal the square root spreading feature for 



E. The Photon Number Distribution 

A complication of random and quantum walks over dif- 
ferent circles is that the number distribution can vary as 
the walker is effectively moving nearer and farther from 
the origin in phase space with the application of each 
generalized Hadamard transformation (|3.ip . This hop- 
ping is responsible for the narrowing of individual peaks 
in Fig. [5] as discussed earlier. 



For 



P(n; N, a, A, A9) = | \p w (N, a, A, A9) \n) \ 

= ^pii^p^NXnlatiN, a, A, A9))(a 3 - {N, a, A, A9)\n) + q t (N)q*(N)(n\f3 t (N, a, A, A9))(p j (N, a, A, A0)\n) 

(3.10) 

the number distribution, the walker's effective distance from the origin in phase space is given by 



][>P(n), 

\ n=0 



(3.11) 



where N, a, A, and AO are suppressed from the expression for brevity, and the walker's radial spread in phase space 
is given by 



Sn=y(h-n). (3.12) 

The expression for n is 

n(N, a, A, A6) = J2 ^ Wp? ( n M ( N > a > A > Ae ) a i ( N , «> A > M ) + * ( N )lj ( N )fi ( N > <*> A > Ad )fo ( N > A > A0 ) > ( 3 - 13 ) 



1,3 



which can be approximated by 

n(N, a, A, AO) w --{ - (a + A) 2 + a(a + A)cosA6> + X(a + A) 
for large N. 



cos NA9 



tAcos [A0(1-7V)] 
sin 2 (A9/2) 



}. (3.14) 



Eqs. (f3TT3|) and ((3TT4|) for the mean number and (|3~T2|) 

for the spread of the walker quantify the degree of hop- 
ping between circles in phase space, and these expres- 
sions will be useful in the next section. Although there is 
hopping to different circles, the QW is clearly evident in 
the quadratic enhancement of phase spreading, with re- 
spect to the Holevo standard deviation, shown in Fig. [3] 



Thus, provided that the parameters a, A, and AO are 
chosen judiciously, the generalization of the Hadamard 
coin flip transformation from (|2.2[) to (|3.1[) does not de- 
stroy the QW, but it docs modify the QW from being 
on a circle in phase space to being on circles in phase 
space. In the next section, we approach the generalized 
Hadamard transformation from the microscopic perspec- 
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tive, and the mean number n turns out to be important 
with respect to controlling the QW in order to ensure 
optimal enhancement of phase spreading. 



IV. FROM JAYNES-CUMMINGS EVOLUTION 
TO QUANTUM WALKS 

In the previous section, we treated the indirectly driven 
coin via the generalized Hadamard transformation (|3.1|) . 
but this transformation was introduced by fiat. In this 
section we consider the JC model Hamiltonian [l2| , which 
underpins so much of quantum optics and cavity quan- 
tum electrodynamics, as a foundation for obtaining the 
generalized Hadamard transformation, or at least a good 
approximation to this transformation under reasonable 
conditions. 

In quantum optics, the simple harmonic oscillator is 
typically the single mode electromagnetic field within the 
cavity, and the coin is an atom transiting the cavity. Cav- 
ity quantum electrodynamic realizations o f Q Ws on the 
circle in phase space have been suggested [ljj EH • 



A. Driven Jaynes-Cummings Model With Large 
Detuning 

For a simple harmonic oscillator with angular resonant 
frequency ui T , coupled with strength g to a qubit of angu- 
lar resonant frequency w a , the JC dynamics for the joint 
system is given by [12( 

Hj C = to? (n + 1/2) + ya a + g(tf&- + aa+). (4.1) 

The joint system is driven by a time-dependent driving 
force (or field) by directly driving the simple harmonic 
oscillator according to 



H dl = e{t) (tie 



t„-io) d t 



ae 



(4.2) 



with e(t) the amplitude and u d the driving carrier fre- 
quency. For simplicity we let e(t) be a constant e for 
some of the time and zero for other times. 

For large detuning g <C |A| = |oj a — u> r \, conjugating 
the JC Hamiltonian under the action of 



(&^ct_ — a<r+) 



(4.3) 



\ ' = cxp 
yields the effective Hamiltonian 

~{u ll + X ^z)n + ^{cJ a + x)^z+0(x 2 ) (4.4) 
for x = g 2 1 'A; the conjugated driving Hamiltonian is thus 



Hdr = VH dr V^ « c (t) 



a + -a. 



he 



(4.5) 



for 'he' designating the Hermitian conjugate. The time 
evolution of Eq. (|4.5[) leads to the generalized Hadamard 
transformation (|3.ip . 



B. Implementation of The Generalized Hadamard 
Transformation 



To implement a QW, first we turn on the driving force 
(e(t) = e) for the Hadamard transformation. In a frame 
rotating at the drive frequency u> d , the effective Hamil- 
tonian of the coin+walker system is thus 

1 



H ct f =- [2x (n + 1/2) - <5da] v z - S dl n 



(4.6) 



with the detuning of the qubit transition frequency from 
the driving force 

<5da = ^d-w a , (4.7) 
the detuning of the resonator from the driving force 

<5 dr = cj d - w r , (4.8) 
and the Rabi frequency 

n R = 2ge/S dl . (4.9) 



The first term in Eq. (J476J) expression effects the coin- 
induced walker phase shift. The unitary operator gener- 
ated by the effective Hamiltonian H e ff is 



exp -i-ffef^H = (if <8> D)E 



(4.10) 



which is a good approximation to the generalized 
Hadamard transformation in (|3.1[) for D(a = — iety) the 
displacement operator (|1 .3[) and 5 is a 'small' operator 
explicitly shown in Eq. (|4. 1 5|) . 

Choosing [2(| co d = 2nx + u; a , H c s then generates rota- 
tions of the qubit about the x axis with Rabi frequency 
Or. In particular, choosing 



and 



co d = 2nx — 2ge/A + u> a 



tn = 7r/20 r 



(4.11) 



= — [A + 2n X -2ge/A] (4.12) 
generates the Hadamard transformation for the coin state 

H = c i *HO R /2* x ^ 13 j 

within the generalized Hadamard transformation (14.101) . 
The choice of pulse duration tu is critical in effecting a 
Hadamard transformation, but this duration itself is a 
function of h, which we know from the previous section 
is time-dependent because the walker is hopping between 
circles in phase space. Specifically in depends inversely 
on Or (|4.9[) . which is itself inversely proportional to the 
driving field detuning 8 dT (|4.8p . The driving field detun- 
ing is a function of u) d (|4.11[) . and uj d is dependent on n 
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(|4.1ip . Therefore, the duration of each pulse in must be 
chosen in accordance with the value of n for the system. 

In order to choose the appropriate pulse duration tu 
for each step, we employ the following protocol, which 
depends on the time-dependent mean number n. In this 
protocol, n is obtained from a theoretical analysis rather 
than continuous measurements or sampling, which could 
disturb the system. In the first step we let 

n=\a\ 2 (4.14) 

and use this value to determine cu^ according to 



Eq. (|4.1ip . Then this value of ui^ is used to compute 
Or and, from this, The duration of the general- 

ized Hadamard pulse is precisely this value of in- In 
subsequent steps n will have changed due to the walker 
hopping to other circles in phase space, so n has to be 
computed and used in a protocol described in Subsec- 
tion IV.D. 



We now have expressions for H and D in Eq. (|4.10l) 
and require 



n ex p 

n=0 

x exp 



(4.15) 



We can see that 5 is close to unity for our choice 
of parameters; thus Eq. (|4. 10|) tends to the generalized 
Hadamard of Eq. (|3.ip . The spectrum of the operators a z 
and a x is 0, 1, and the relative size of fit + a and i(cv — a) 
is never much more than \a\ because 1(d) | = \a\. In the 
case t-ftX = irg/4e C 1, we neglect the higher orders of 
tnX- In the case of large detuning, that is g/Sdi- <C 1, the 
term 



4<5 d 



(4.16) 



can also be neglected. Thus 3 can be approximated by 



c w exp 



(4.17) 



The evolution of initial states under S are shown as 

H|j,a)w \j,a+{-iy7d K g/8),j = 0,l. (4.18) 



As E in Eq. (|4.17[) is close to an identity operation for 
the restricted choices of parameters, the resultant gener- 
alized Hadamard transformation (|4.10[) is quite close to 
the ideal (|3.ip in the previous section. It is thus impor- 
tant to choose parameters for which 2 can be neglected. 
In this case, the displacement operator D in Eq. (|4. 10[) 
is responsible for displacing the walker's distance from 
the origin in phase space by \a\ t— > |a|(l + etn/2). For- 
tunately, even the effects of this induced jump in |a| can 
be minimized by varying the duration of successive gen- 
eralized Hadamard pulses. 



C. Implementation of The First Step 

In the previous subsection, we have seen how the gen- 
eralized Hadamard transformation generated by the JC 



model is very close to the ideal Hadamard transformation 
of Sec. III. The importance of choosing the appropriate 
duration of the generalized Hadamard pulse was noted 
in Subsection IV. B. Each step of the QW corresponds to 
first performing the generalized Hadamard transforma- 
tion and then the conditional phase shift operation given 
by F (|2.3p . In this subsection we concentrate solely on 
the walker's first step, which is the generalized Hadamard 
transformation followed by F. 

The conditional phase shift A9 has a size that is con- 
strained by (|2.4p . In terms of parameters in the JC 
model, the step size is 



A9 = ± X (T + t H ), 



(4.19) 



for t the time between generalized Hadamard pulses. Be- 
cause the JC Hamiltonian applies to the dynamics both 
during the generalized Hadamard pulse, which has du- 
ration tfj, and during the period between these pulses, 
which has duration r, the step size (|4.19p is proportional 
to the total time for each step, namely r + fn- 

At time r + in the first step is completed, but n has 
changed. The new n after the completion of the first step 
is required to calculate the appropriate tn for the second 
step. The value of n after the first step is readily obtained 
from Eq. (|3 . 13[) by inserting the relevant parameters as 
well as N = 1. From this value of n, the pulse dura- 
tion for the next generalized Hadamard transformation 
is given by Eq. (|3 . 13|) . This knowledge of £h for the next 
generalized Hadamard transformation prepares us for the 
second step. 
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D. N Steps 

The previous subsection describes how to perform the 
first step and obtain the information required to set the 
duration for the subsequent generalized Hadamard trans- 
formation. In this subsection we describe the transforma- 
tions required for the walker to go an arbitrary number TV 
steps. Unlike the case of the quantum walk on a single 
circle or the case of quantum walks on circles described 
in Sec. Ill, here the choice of £h for each circle is more 
complicated but quite important. 




FIG. 4: (a) The Holevo Standard deviation of phase for both 
the quantum and classical random walks up to TV = 15 for a = 
3, d = 21, and typical system parameters (w a , ui r , g, e)/2ir = 
(7000,5000,100,1000) MHz [23]. Numerical simulations re- 
veal that the Holevo standard deviation is almost independent 
of the initial state of charge qubit, and is approximately linear 
in TV, an = (1.3964 ± 0.0180)t + (0.1208 ± 0.0018). (b) The 
numerically simulated Holevo standard deviation for phase 
distribution in log-log scale is shown to be approximately lin- 
ear in logt: logo-H = (0.924 ± 0.009) logt + (0.442 ± 0.004) 
and the r coefficient is 0.99. 

For an arbitrary i th step, we can calculate the average 
photon number n(i, a, A, A9) based on the analytical re- 
sult (|3.13[) . and then decide the pulse duration of the i th 
step (|4~T2|) 



Age 



[A + 2n(i,a,\,A6)x-2ge/A]. (4.20) 



We apply the generalized Hadamard transformation 



cxp 



d-ffcff^H 



following by the unitary operator of the 

free evolution exp \}x( T + *h)1 n&z- These two applica- 
tions together effect unitary operation 



U eS w F(H ® D) 



(4.21) 



Using our protocol for choosing durations of general- 
ized Hadamard pulses, we obtain numerically the Holevo 
standard deviation for the phase distribution of the re- 
duced walker state as a function of time t. In contrast 
to the related plots in Fig. [3] of Sec. Ill, which depend 
on the number of steps TV, these plots explicitly depend 
on t. In Sec. Ill the choice of TV vs t is not significant 
because t cx TV; here, however, t is not proportional to TV 
because of the varying duration of each step due to the 
variability of n. In physical systems, the random walk is 
characterized by its time dependence so, in that spirit, we 
also use time t, rather than the number of pulses TV, to 



show the quadratic enhancement of the phase spreading 
for the QW vs the random. 

This quadratic enhancement is evident in Fig. |4] To 
show this more explicitly we apply linear regression tech- 
niques to the log-log plot, which theoretically should be 
linear with a slope of 1/2 in the classical case (depicted 
as a solid line in Fig. [4][b)) and slope 1 in the quantum 
case for small spreading of phase. The linear regression 
results are presented in detail in the caption of Fig. [4] 
and residual r = 0.99 for the QW, confirming the linear 
relationship between cth vs t. The slope is 0.924, which 
is quite close to unity. Together the slope being close 
to unity and the high value of r demonstrate that this 
protocol does indeed lead to an enhancement of phase 
spreading that is very close to quadratic and is thus a 
signature of QW behavior. 



V. CONCLUSIONS 

Motivated by the physical difficulty of directly driv- 
ing a coin qubit in a cavity quantum electrodynamical 
realization of the quantum quincunx, we generalized the 
Hadamard coin flip to also kick the resonator. In Sec. Ill 
this kick was incorporated within an idealized generalized 
Hadamard transformation, and Sec. IV approximately 
obtained the generalized Hadamard transformation di- 
rectly from the ubiquitous Jaynes-Cummings Model. 

The generalization of the Hadamard transformation 
modifies the walk from being on one circle in phase space 
to hopping between circles in phase space. Despite this 
hopping, the quantum walk is evident, in the quadrat- 
ically enhanced spreading of phase. In Sec. IV the du- 
ration of each generalized Hadamard pulse is modified 
according to which circle the walker is on — equivalcntly 
the time-dependent mean number n — which means that 
the spreading of phase in time t is slightly different from 
spreading as a function of number of steps TV. We show 
the quadratic enhancement in terms of the more experi- 
mentally relevant time t, which is the signature of quan- 
tum walk behavior. 

As explained in (l3j . the quantum walk behavior can 
be ascertained by bringing in controllable decoherence. 
Then tuning of decoherence will interpolated the phase 
spreading from linear in time to the square root of time. 
Furthermore, although phase is not directly measured, 
its cosine and sine can be inferred from homodyne mea- 
surements, or from full optical homodyne tomography. 

Appendix 

We calculate the state of the coin+walker system after N 
steps for N small. For the initial state |<&) in Eq. (|3.5p . 
after the TV th step of walking on the circles, the state 
|$(TV)) = U N |$) is shown in Eq. (|3T6]l . where the coef- 
ficients ft (TV) and qi(N) are obtained from the following 
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recursion relations (for N > 2) For the case N = 0, we have Po(0) — and <7o(0) 

— i/v2- We will show the case iV = 1 below. 

f K (7V-l)/V2 ifl<z<2^ 2 
M j " I ft- 2 ^-»(JV - 1)/V2 if 2 W - 2 < i < 2 N ~\ 



(Al) 



and 



% (JV) = \ Pi ^ N ^ 1 n-1~ 2 n-1 The coherent state with a* (A^ A, A6>) and &(N, A, A6>) 

| — Qi-2 N - 2 (N — l)/v2 if 2 N < i < 2 N 1 , carL gj so obtained from the following recursion relations 

(A2) for N > 1 



a^A^Atf) 



N-2 



[ati(N- l,A,A(9) + A]e iAe if 1 < i < 2 

[/3 l _ 2 «- 2 (A^ - 1, A, A9) + A] e iAe if 2 W - 2 < * < 2 N ~\ 



(A3) 



and 



[a>i(N- l,A,A6») + A]e- iAe 
[A-2^(iV-l,A,A0)+A]e- 



if 1 < < < 2 W ~ 2 
if 2 W ~ 2 < » < 2 X 



(A4) 



For the case N = 0, a (0) = /3 (0) = a. 

After the first step, the state of the system is 

|*(1)) =p 1 (l)|0,ai(l)> + «i(l)|l l) 9i(l)) 1 (A5) 

with 

Pi(l) = -^,gi(l) = — , (A6) 
ai(l) = (a + A)e iAe , /3r(l) = (a + A) e - iAe . 
After the second step, the state is 

2 

| $(2)) = J>(2) |0, a, (2)) + gi(2) |1, A(2)> , (A7) 

i=l 

with 

" w -^" w -^ (A8) 

1 + i , - 1 — i 
^) = ^f.*(2)— ^, 

ai (2)=ae 2iA9 + X(e 2iAe + e iAe ), 
a 2 (2) = a + A(e iAe + 1), 
ft (2) = a + A(l + e- iAe ), 

o /ol -2iA9 i \ / -iA0 i -2iA0\ 

p 2 (2) = ae + A(e +e ). 
The third step leads the state to 

4 

|*(3)>=£ft(3)|0,a i (3)> + g i (3)|l,A(3)>, ( A 9) 



with 



Pi(3) = 
Ps(3) = 
9i (3) = 



1 


+ i 




4 


1 


+ i 




4 


1 


+ i 



,Pa(3) 



1 - i 



,P4(3) = - 



4 ' 
1 - i 



(A10) 



92(3) 



4 
1 - i 



4 4 . 

9s(3) = 4-^*1(3) = - 4 -^, 

a! (3) = ae 3iAe + A(e 3iAe + e 2iAe + e iAe ), 
a 2 (3) = ae iAe + A(e 2iAe + 2e iAe ), 
03(3) = ae iAe + X(2e iAe + 1), 



a 4 (3) 


= ae 


-iA0 


f A(e iA0 + l + e- iAe ) 


Pi (3) 


= ae 1 


Ae , 


A(2e iAe + e" iA9 ), 


/3 2 (3) 


= ae 


-iA0 


f A(l + 2e- iA9 ), 


ft (3) 


= ae 


-iA0 


f A(2e- iA9 + e~ 2iAe ), 


ft (3) 


= ae" 


-3iA0 


+ A(e" iAe + e" 2iA9 4 



-3iA0\ 



The entanglement between the coin qubit and the su- 
perposition of coherent states leads to the signature of 
QW compared to random walk, that is the quadratic in 
phase spreading. From Fig. [5l for the given a and fixed 
A9, the hopping between circles, i.e. A leads the phase 
distribution to be skewed towards positive <fi and indi- 
vidual peaks can become narrower or broader. However, 
for the case A <C a, we still obtain the characteristic 
quadratic enhancement in phase spreading for QW. 
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FIG. 5: The phase distribution for the walkers' location after 
the first three steps of the QW on the different circles with 
initial state (|0) + i |1)) \a = 3) /y/2, A6 = 0.35 and different 
A. (a) A = and (b) A = 0.4. The yellow line is for the case 
N=0, the blue one for N=l, the green one for N=2 and the 
red one for N=3. 
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